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Summary 

The objective of this investigation was to validate 
two propagation models (modal decomposition and 
finite element) for predicting the performance at low 
frequencies of duct liner configurations built up from 
an array of seven individually tuned resonators. Most 
tests were conducted on a configuration consisting 
of round resonator cavities. However, limited tests 
were also conducted on a configuration with square 
resonator cavities. These analytical models had pre- 
viously been developed for high frequency attenua- 
tion in aero-engine inlets, but no satisfactory valida- 
tion of their accuracy existed in the low frequency 
range. In the low frequency application, the liner 
surface impedance is characterized by spatial varia- 
tions as a consequence of available structural compo- 
nents and fabrication techniques. The scale of these 
variations may cause the liner to appear acoustically 
nonuniform. 

The modal decomposition model was of the 
“infinite waveguide” type and consequently assumed 
a uniform surface impedance distribution. The 
finite element model, on the other hand, could handle 
surface impedance nonuniformities in the axial direc- 
tion but assumed uniformity in the transverse direc- 
tion. Thus, to provide acceptable impedance bound- 
ary conditions for the two predictive models, two pro- 
cedures were used to modify the impedance distribu- 
tions. A “smeared” impedance distribution removed 
spatial nonuniformities by equating volume veloci- 
ties between the resonator active area and its respec- 
tive apportioned area in the array. This procedure 
resulted in a specific impedance change. An “equiv- 
alent length” impedance distribution procedure sim- 
ply distorted the active circular area of each res- 
onator into an equivalent length, such that the 
specific impedance was preserved. For the modal de- 
composition model, both procedures were applied to 
provide a uniform surface impedance free of trans- 
verse and axial dependence. The equivalent length 
distribution procedure was used to remove only 
the transverse impedance variation so that a two- 
dimensional finite element propagation model could 
be implemented with axial impedance variations left 
intact. 

Attenuation rates predicted with the modal de- 
composition propagation model were generally in 
good agreement with measured attenuation rates, for 
the round and square resonator arrays. This was 
true for both impedance distribution modeling pro- 
cedures and for frequencies up through the first reso- 
nance of the resonator backing cavities, or for nondi- 
mensional frequency parameter kh values of 0.8. 
Attenuation rates obtained from the finite element 


model were in excellent agreement with the measured 
values, as were the attenuation profiles. For the 
square resonator array, the predicted and mea- 
sured attenuation profiles did not agree as well, al- 
though the attenuation rates were still in very good 
agreement. 

Predicted and measured transmission losses for 
test frequencies beyond the resonator cavity first res- 
onance were in very good agreement when all res- 
onator cavity depths were set the same. For linearly 
varying cavity depths, predicted and measured trans- 
mission loss trends were still in good agreement, but 
less so than for uniform cavity depths. 

Detailed calculations of attenuation and phase 
profiles, at axial locations near the resonator array 
elements and in the presence of mean flow, revealed 
a fine-scale structure of acoustic “hot spots” superim- 
posed on the general acoustic field structure. Care- 
ful measurements verified the existence of such hot 
spots, which apparently are attributable to convec- 
tive effects alone, since no flow-induced impedance 
changes were assumed in the model. 


Symbols 

A 

amplitude coefficient 

Amn i B mn 

amplitude coefficients 
associated with the mnth 


acoustic mode 

C 

sound speed 

dB h 

attenuation per duct height, 
in dB 

D 

single resonator diameter 

h 

duct height 

i 

= x/— T 

I 

acoustic intensity 

h 

incident acoustic intensity 

It 

transmitted acoustic 
intensity 

[J] 

Jacobian matrix 

k 

free-space acoustic wave 
number, u/c 

kh 

nondimensional frequency 
parameter or nondimen- 
sional duct size 

kmn 

total wave number associ- 
ated with the mnth acoustic 
mode 

l 

typical duct dimension 

L 

single resonator length 


L e 

axial dimension 

m. n 

spinning wave numbers 

M 

centerline Mach number 

N\. A ’2 . A, 3 . N4 

finite element shape 
functions 

p(x.y.z) 

acoustic pressure 

PA 

acoustic pressure at leading 
edge 

PB 

acoustic pressure at trailing 
edge 

Ps 

static pressure 

Pt 

total pressure 

P12 

pressure ratio, 

Q 

constant (see eq. (7)) 

r 

duct radius 

ri,Xi 

global coordinates of corner 
nodes 

R 

reflection coefficient 

Re 

Reynolds number 

R f 

flow resistance 

t 

time, e~ lojt time convention 
assumed 

T 

temperature, °C 

U, V, w 

acoustic particle velocities 
in x-, y-, and ^-directions, 
respectively 

Ui,U2,U3,U4 

values of u , v, w, or p at 
the four node points at the 
numbered corners of the 
element 

V* 

complex conjugate of 
acoustic velocity 

W 

resonator width 

X 

distance measured along 
duct axis 

xi, X2 

axial positions 

y 

height above lower wall 
of duct (resonator array 
surface) 

Ym(y )i z n (z) 

eigenfunctions 


z 

distance across duct in 
horizontal direction 

Z 

resonator acoustic 
impedance 

z E 

equivalent acoustic 
impedance 

Zs 

smeared acoustic impedance 


wave number in 2 -direction 

Pi, fa, fa P 4 

wall admittances 

1m 

wave number in y-direction 

f 

acoustic impedance, normal- 
ized by pc 

v,Z 

set of local coordinates 
within the region 

p 

density 

p 

mean density 

OJ 

source frequency 

V 

Laplace operator 

Abbreviations: 


LE 

leading edge 

SPL 

sound pressure level 

TE 

trailing edge 

Introduction 



The efficient suppression of unwanted sound in 
ducted air-moving systems is of continuing interest 
for numerous applications. Although silencer tech- 
nology is very advanced for many applications (such 
as automotive exhaust mufflers and turbojet engine 
inlets), new and unique requirements continue to 
make severe demands on current technology. One 
such requirement arises from the urgent need for 
ultraquiet wind tunnels for specialized aerodynamic 
tests, such as the investigation of transition from lam- 
inar to turbulent flow in the absence of acoustic ex- 
citation. The receptivity of laminar flow to acoustic 
excitation is of particular interest, since it may be 
a barrier to obtaining the significant drag reduction 
associated with laminar flow (see refs. 1, 2, 3, and 4). 

Abrahamson (ref. 5) proposed an ultraquiet wind- 
tunnel design for laminar flow transition research 
which set an acoustic noise level limit not to ex- 
ceed 60 dB in any octave band from 31 to 4000 Hz. 
Such a requirement is particularly difficult to meet at 
low frequencies because of turbulent processes such 
as aerodynamic throttling, interaction with vibrat- 
ing panels, and flow turning. One way to obtain 
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greater attenuations at low frequencies is to inte- 
grate individually tuned resonator absorbers into a 
“smooth” wall diffuser. Thus, the sound absorbing 
walls must be aerodynamically smooth and conform 
to the desired streamline flow pattern to minimize 
self- noise and to maintain flow quality. Tradi- 
tional automotive- type, low frequency resonator muf- 
fler technology, which incorporates flow direction 
changes and flow path discontinuities to attenuate 
exhaust noise, is unacceptable for these stringent 
requirements. 

The development of resonant liner technology was 
greatly stimulated in the early 1970’s by the intro- 
duction of the turbofan engine into the civilian air- 
line fleet. An intense research effort was directed to- 
ward reducing high frequency noise (2.5 kHz to about 
10 kHz) radiating through turbofan inlet ducts to the 
environment. Partly because of size, weight, and per- 
formance advantages, resonant duct liners with rela- 
tively smooth surfaces were the predominant choice 
for duct liner acoustic treatments. Lowson, in ref- 
erence 6, gives an overview of duct liner technology 
with emphasis on application to aero-engine ducts. 
Similarities between the geometry of the aero-engine 
duct and that of the present application suggest 
the applicability of the same technology, but at 
wavelength scales 10 to 100 times larger. 

The purpose of this paper is to describe a 
scale model, experimental validation of two theo- 
retical techniques available for use in low-frequency- 
attenuating duct liner applications. One aspect to 
be investigated was the surface impedance distribu- 
tion modeling procedure for specifying the bound- 
ary conditions to achieve acoustic attenuation rates 
predicted by these techniques. To this end, two 
different scale model resonator arrays were fabri- 
cated, and their acoustic performances were mea- 
sured and compared with analytical predictions us- 
ing both a closed form modal approach and a finite 
element analysis. The following sections of this re- 
port describe comparisons between the two theoret- 
ical models on the basis of calculated and measured 
axial attenuation rates, attenuation and phase pro- 
files, and transmissions losses for a scale model, lin- 
ear array of resonators. Individual resonators con- 
sisted of round backing cavities fitted with mov- 
able pistons for variable tuning. The resistive 
element was provided by a porous fiber metal screen 
bonded to the resonator openings. The resulting 
resonator array was thus characterized by circular 
regions of relatively low impedance separated by 
regions of “hard wall.” Limited tests were also con- 
ducted on a resonator array with square backing cav- 
ity cross sections, in which case the hard wall inter- 
stitial regions between resonators were considerably 


reduced compared with those for the round resonator 
array. 

Experimental Setup and Measurement 
Procedures 

Test Models 

Figures 1 and 2 show sketches of the round and 
square cavity resonator array test models, respec- 
tively. The structure of the round cavity resonator 
array of figure 1 will be discussed in detail. The 
same remarks also apply to the square cavity res- 
onator array, except for the obvious geometrical dif- 
ference. The seven resonator cavities were fabricated 
from 1.75- in-id aluminum tubing with a wall thick- 
ness of 0.125 in. To form a linear array of seven 
resonators, the cylindrical cavities (round or square) 
were set into an epoxy fixture as shown in the figure. 
Cavity depths were adjustable from 0 to 5.4 in. by 
means of adjustable pistons equipped with O-ring 
seals. It should be noted from the figures that the 
overall axial length and width of both resonator ar- 
rays were 15.5 in. and 2.0 in., respectively. Thus, 
each resonator could be apportioned an axial length 
of 2.2 in. The resistive component of the resonator 
impedance was supplied by a porous fiber metal face 
sheet of thickness 0.028 in., with a nominally rated 
flow resistance of 100 kg/m 2 -sec. A single strip of this 
fiber metal, with dimensions 15.5 in. x 2.0 in., was re- 
moved from a parent sheet and bonded to the epoxy 
fixture. Thus, the resistive elements for each res- 
onator consisted of circular regions of the face sheet 
congruent with the resonator mouth. 

Figures 3(a) and (b) show photographs of the two 
completed resonator arrays. Figure 3(a) shows a view 
of the arrays from the front, or active, side, and fig- 
ure 3(b) shows a view of the backing cavities with a 
few of the airtight, adjustable pistons inserted. The 
square cavity resonator array is shown attached to 
the fixture that allows the resonator array to be in- 
stalled on the bottom wall of the flow impedance 
tube apparatus. Note that in figure 3(a) there ap- 
pears to be some “wicking” of the bonding agent into 
the resonator mouth areas. This effect appears most 
pronounced for the square resonator array. How- 
ever, visual inspection, by observing differences in 
light transmission through the face sheet over the 
resonator mouth areas, indicated that the square res- 
onator active face-sheet areas were evidently more 
uniform than those of the round resonator array. 

Resonator Normal Incidence Impedance 
Measurements 

Although the contribution of the resonator cavity 
reactive component to the total impedance can be 
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reliably calculated, such is not the case for the re- 
sistive face-sheet contribution. The manufacturer’s 
nominal flow resistance value is merely a selection 
guide and may be subject to a variance of up to 
25 percent from sample to sample. Thus, to ascertain 
a more accurate estimate of the acoustic impedance, 
the normal incidence impedance was measured for a 
typical resonator from each resonator array and flow 
resistance measurements were conducted on all the 
resonator resistive face sheets. 

Normal incidence impedance data were obtained 
by means of the two-microphone method, imple- 
mented with one of two “impedance tubes” (one of 
which is shown in fig. 4). These impedance tubes 
were constructed of 1-ft lengths of the same cylin- 
drical tubing used for the resonator array cavities, 
thereby avoiding area discontinuities between the 
impedance tube and the active resonator area. Han- 
dles were attached to hold the device in place against 
the resonator whose impedance was to be measured. 
A seal between the tube and the resonator face 
was provided by a rubber gasket, as shown in the 
photograph. 

Instrumentation for the two-microphone imped- 
ance measurement method is shown in figure 5. 
Acoustic plane waves in the frequency range 0.5 
to 3.0 kHz were propagated down the tube and 
were normally incident on one of the seven res- 
onator array elements, which was aligned and pressed 
against the gasketed, open end of the impedance 
tube. The amplitude and phase for each frequency 
of interest were measured at two probe ports lo- 
cated 2.08 in. and 3.08 in. from the face of the test 
specimen. The microphones used in this measure- 
ment were ^-in. condenser microphones connected 
to j^r-in-id stainless steel probe tubes, with lengths 
such that the ends were flush with the impedance 
tube inside wall. Amplitude and phase data were 
acquired by way of the signal path illustrated at 
the left of the diagram. This data acquisition sys- 
tem was essentially the same as that used to acquire 
the axial attenuation and phase profiles and will be 
described in more detail in the “Acoustic Measure- 
ments in Flow Impedance Tube” section. To en- 
sure constant levels at the resonator face during 
impedance measurements, signals from a |-in. refer- 
ence microphone, located adjacent to the resonator 
face, were input to a feedback loop (automatic level 
regulator). 

Resonator Flow Resistance Measurements 

The dc flow resistance of a thin porous mate- 
rial is defined as the pressure drop across the ma- 
terial per unit incident steady flow velocity. It is 
a useful indicator of the acoustic resistance if the 


test specimen is thin compared with the relevant 
acoustic wavelengths and if flow velocities are not 
too high compared with the acoustic particle veloci- 
ties of interest. To conveniently measure and evalu- 
ate the variability of the acoustic resistance values of 
all the resonators, a computer-controlled valve and 
flowmeter arrangement was used to measure the flow 
resistance of each resonator. To accomplish this mea- 
surement, steady airflow was forced through the face 
of each resonator, by way of the backing cavities, at 
a series of precisely controlled velocities. The princi- 
ple of operation of such an apparatus is described in 
reference 7. Corresponding pressure drops and flow 
velocities across the resistive face sheet of each res- 
onator were measured and used to compute the flow 
resistance. Clearly, flow resistance is far more conve- 
nient to measure than acoustic impedance. Thus, in 
this investigation, flow resistance measurements were 
used to determine the acoustic resistance variations 
among the resonators. 

Flow Impedance Tube 

In this experiment, several different acoustic tests 
were performed on the resonator arrays to ascertain 
their acoustic characteristics in a grazing incidence 
and/or grazing flow environment. The tests were 
conducted in the flow impedance test apparatus at 
the Langley Research Center, as shown in figure 6 
(ref. 7). This multiconfigurational apparatus is de- 
signed to produce a controlled aeroacoustic environ- 
ment with sound pressure levels up to 140 dB and a 
flow speed of up to Mach 0.6, over a test specimen 
length up to 15.5 in. Moving from left to right in the 
figure, it is convenient to divide the apparatus into 
three sections as indicated. The source section begins 
with a supply plenum from which air (if flow is de- 
sired) percolates into a 2.25-in-diameter round duct 
constructed of high-resistance fiber metal. A ther- 
mocouple located at the plenum exit provides a con- 
trol point for airflow temperature, and a pitot-static 
tube located about 12 duct diameters downstream 
from the plenum exit provides a means for determin- 
ing a reference centerline Mach number. An electro- 
magnetic, 120-W acoustic driver is used to propagate 
either broadband or pure tone plane waves through 
the square cross section (2 in. x 2 in.) of the test 
section region. The acoustic waves are gradually 
combined with the flow through the high-resistance 
(60000 kg/m 2 -sec) cylindrical duct as shown. Cross- 
section areas are maintained constant throughout the 
system. Also, over a region from 5 in. upstream 
of the test linear leading edge (LE) to 2 in. down- 
stream of the test linear trailing edge (TE), axially 
uniform flow conditions can be maintained. 
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Over the test specimen region (indicated by LE 
and TE in the figure), and 48 in. upstream of it, a 
0.5-in.-wide precision-machined slot in the top wall 
allows an axial bar to traverse the test section length 
by means of a digital stepping motor under computer 
control. The traversing bar will accommodate the 
installation of several different types of transducers. 
In this experiment, measurements were made with 
either a flush-mounted condenser microphone or a 
strain gauge pressure sensor incorporated into a spe- 
cially designed, streamlined probe, as described in 
reference 8. Under control of data acquisition pro- 
grams, these sensors could sample the acoustic field 
at desired locations. These data provided the at- 
tenuation and phase profiles over the test specimen 
length (as suggested by the sketch beneath the test 
specimen location) and information for calculating 
acoustic intensities upstream of the leading edge and 
downstream of the trailing edge. 

Test specimen lengths up to 15.5 in. can be ac- 
commodated in the bottom wall of the test section. 
With the full length of the test specimen installed, a 
transducer mounted on the traversing bar can probe 
the entire length of the test specimen, plus about 
48 in. upstream of the test specimen leading edge and 
at least 2 in. downstream of the trailing edge. In this 
manner, the acoustic wave fields upstream and down- 
stream of the test specimen can also be measured. 
A more detailed description of this axial-traversing 
measurement system is given in reference 8. For 
some tests, total attenuations across the resonator ar- 
ray were obtained with |-in. condenser microphones 
flush-mounted in the side walls of the tube at the 
leading and trailing edges. 

Acoustic waves exiting the test section region are 
absorbed by an anechoic termination section spe- 
cially designed to be nonreflective with or with- 
out mean flow present. This section consists of a 
20-ft-long 2-in. x 2-in. duct, with the two absorptive 
vertical walls consisting of a 30-percent perforated 
plate backed by open-celled foam compressed to one- 
seventh the original volume at the entrance end of the 
termination. The compression gradually decreases 
over the 20-ft termination length. A 2-ft foam wedge 
is encountered near the termination end, as shown 
in figure 6. This design was intended to inhibit flow 
noise and to allow the acoustic waves to be absorbed 
without significant reflection back into the test 
section region. 

Airflow through the flow impedance tube was pro- 
vided by the combination of a pressurized plenum 
in the source section and a vacuum source at the 
exit of the termination section. An automatic 
electromechanical control system maintained con- 
stant flow through the test section to within about 


3 percent of a preset Mach number. Static pres- 
sure at the leading edge of the test section was 
maintained to within about ±3 mm Hg of a preset 
value, and the free-stream static temperature of the 
airflow was maintained to within ±2°C of a preset 
value by means of electric heaters in the air supply 
line. From experience, flow stability was optimum 
when the static pressure at the test section was main- 
tained 5 to 10 mm Hg below local barometric pres- 
sure and the static temperature was maintained at 
ambient temperature. 

Acoustic Measurements in Flow Impedance 
Tube 

Acoustic tests conducted in the flow impedance 
tube apparatus consisted of broadband attenuation 
measurements, axial attenuation and phase profile 
measurements for pure tones, and transmission loss 
measurements. The instrumentation for accomplish- 
ing these measurements is depicted in figure 7. De- 
pending on whether a pure tone or broadband exci- 
tation was desired, the power amplifier was driven by 
a programmable frequency synthesizer or a random 
noise generator. To obtain the broadband attenu- 
ation measurements across the resonator array, the 
transfer function magnitude was measured between 
the microphones located at the leading and trailing 
edges by means of the spectrum analyzer shown at 
the right of the diagram. To obtain the pure tone ax- 
ial attenuation and phase profiles, the signal process- 
ing path for the presssure probe signal shown at the 
center of the diagram was used. After appropriate 
signal conditioning, the signal was filtered through a 
50-Hz bandwidth tracking filter to remove harmonics 
of the signal generated by the acoustic driver, which 
was operated near its peak capacity. To further im- 
prove the dynamic range of the measurement process, 
the signal was then subjected to synchronous averag- 
ing to remove residual random noise. This step was 
really only necesary with flow, but was also included 
in the no-flow measurements for consistency. Abso- 
lute sound pressure level and phase were then calcu- 
lated with the computer by using stored calibration 
information. These calculated values were also stored 
for later use to generate axial attenuation and phase 
plots. 

Transmission losses were obtained by determining 
acoustic intensities / t - and It associated with the in- 
cident and transmitted plane waves in the hard-wall 
regions of the test section upstream and downstream 
of the test specimen leading and trailing edges. At 
the leading and trailing edges of the resonator array, 
the plane wave is necessarily contaminated by higher 
order modes to provide continuity of the acoustic 
field across the discontinuity between the hard and 
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soft walls. These higher order modes are below their 
cut-on frequencies, however, and tend to attenuate 
rapidly with distance from the discontinuity. Thus, 
to obtain valid measurements of transmission loss for 
the plane wave, the acoustic intensities should be ob- 
tained at sufficient distances from the leading and 
trailing edges to ensure that higher order modes have 
decayed. Experimentally, this will be the case if the 
intensity measurements are location independent. 

The intensity measurement procedure, from 
which transmission losses, 101og(/ 2 //f), were ob- 
tained, is described in the “Analysis” section and 
can be treated as a modification of the two- 
microphone impedance measurement method. In the 
present situtation, the computer-controlled micro- 
phone positioning capability allowed a single micro- 
phone, flush mounted on the traversing bar, to re- 
place the two microphones normally used in the two- 
microphone method. In this manner the need for 
accurate, and tedious, microphone calibrations was 
eliminated. Measurements were made 10.0 in. and 

9.0 in. upstream of the liner leading edge, and 

2.0 in. and 1.0 in. downstream of the liner trailing 
edge. The data acquisition was done in a manner 
almost identical to that used in the normal incidence 
impedance measurements. Since the results were 
not position dependent, direct comparisons between 
the predicted and measured transmission losses could 
then be made. 

Analysis 

Requirements and Approximations 

It is well known that acoustic propagation phe- 
nomena through lined ducts scale on the dimension- 
less parameter kh, where k = uj/c and h is a typ- 
ical duct dimension. This fact permits propagation 
model validation experiments performed over one fre- 
quency range to be applicable to another frequency 
range. The test results presented in this paper were 
obtained over a kh range of 0.3 to 2.5 (0.3 to 3.0 kHz), 
where h is the duct height. 

A rigorous, analytical treatment of sound propa- 
gation in a muffler that is also designed to double as 
a diffuser and that incorporates individually tunable 
resonators distributed over the walls, would require a 
fully three-dimensional (3-D) analytical model. The 
importance of 3-D effects, however, depends on many 
factors, including acoustic source frequency, duct ge- 
ometry and its rate of change with axial distance, 
and acoustic wall impedance distribution. Since no 
modeling code was available that would fully incorpo- 
rate all possible 3-D effects, two separate codes were 
used which employed different sets of simplifying as- 
sumptions. The results of both analyses were then 


compared to provide a cross-check on these assump- 
tions. These codes were the modal analysis code de- 
veloped by Watson (refs. 9, 10, and 11), and the finite 
element code ADAM (refs. 12 and 13). A brief de- 
scription of the two methods follows. 

Two-Dimensional Modal Decomposition 
Model 

The disadvantage of the Watson formulation is 
that it is unable to handle the effects of axial vari- 
ations either of duct geometry or of duct liner 
impedance. Thus, it was necessary to employ an 
“equivalent” axially uniform impedance to approx- 
imate the impedance distribution. Two different 
methods of deriving this equivalent impedance were 
used. In an a priori sense, either of these methods 
is equally acceptable, since they both employ an as- 
sumption to smooth out the nonuniformities present 
in the surface impedance distribution. These two 
methods for deriving the equivalent impedance dis- 
tribution will be described with the aid of figure 8. 
Figure 8(a) depicts the actual geometry of the res- 
onator array impedance distribution drawn to scale. 
The circular shaded areas of diameter D represent 
the resonator impedances Z. It should also be noted 
that exactly the same procedure was used for the 
round and square resonator arrays. 

Smeared impedance method. In this modeling pro- 
cedure, each resonator impedance was uniformly dis- 
tributed over its apportioned length L and width W 
in the array, as suggested in figure 8(b). Physically, 
this means that acoustic volume flows through the ac- 
tual and idealized resonator face sheets are equated. 
Thus, the specific acoustic impedance is different for 
the idealized resonator face sheet by the ratio of the 
so-called smearing factor. The smearing factor on 
the measured impedance was 1.83 for the round res- 
onator array and 1.44 for the square resonator ar- 
ray. Both the smeared and equivalent length (to be 
described next) impedance distribution models were 
used with the two-dimensional (2-D) modal decom- 
position propagation model. 

Equivalent length method. With reference to 
figure 8(c), the circular area of a single resonator of 
diameter D in the array is equated to an equivalent 
rectangular area with axial dimension L e and width 
W. Thus, for the circular resonator diameter of 
1.75 in., the resultant equivalent length is 1.2 in. This 
equivalent length corresponds to about 55 percent of 
the apportioned length of 2.2 in. for each resonator 
element and therefore for the entire resonator array 
length. Note that this modeling procedure conserves 
specific acoustic impedance, since it only changes 
the shape of the impedance region. Figure 8(d) 
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shows the equivalent length procedure as applied 
in modeling the array impedance distribution for 
the finite element code to be described later in this 
section. The advantage of this code is the ability to 
account for axial variations in impedance. 

Formulation of two-dimensional modal decom- 
position method. The formulation of the 2-D modal 
decomposition method is based on a closed-form solu- 
tion to the Helmholtz equation in a rectangular duct 
by separation of variables. In the absence of mean 
flow, acoustic waves propagating in a duct (fig. 9(a)) 
are governed by the 3-D Helmholtz equation 

V 2 p + k 2 p = 0 ) (1) 


Therefore, 


Q = 


-ikfi 3 
7m 


and, 


Ym — COS 7rnJ/ 


ikfo . 

sin 7 m i/ 

7 m 


Similarly, at z = 0 for Z n , equation (6) gives 


( 8 ) 


Z n — cos ot n z 


ikfii 

a n 


sin a n z 


(9) 


Substituting these expressions in equations (3) and 
(5) we get the transcendental equations 


A harmonic solution in x is of the form 


OO OO 

p{x,y,z) = ^2 y 2 A mne lkmnX 

m— 1 n = 1 


1 d —ikmnX 

~r E>mn & 

) Y m (y) Z n (z 

0 (2) 

where the eigenfunctions Y m and Z n must satisfy the 
boundary conditions 

^ - ikfaYm = 0 

dy 

' 5 T 

c+ 

II 

( 3 ) 

^ + ikfoYm = 0 
dy 

(at y = 0) 

(4) 

■, „ a ? n 
lkp2^n — 0 
oz 

(at z = l) 

(5) 

+ ikfilZn = 0 

oz 

(at z = 0) 

(6) 

(Wall admittances /? are shown in fig. 9(a).) 
a solution for Y m of the form 

Suppose 

Ym = cos 7 m y + Q sin 7 m y 

( 7 ) 


then 


dY n , 

dy 


= -7 m sin 7 m y + ImQ cos 7 m y 


or, substituting in equation (4) 

ImQ + ikfi 3 


tan 7 m y + 


= 0 


-7 m + ikfoQ 

Now, when y = 0, tan 7 m y = 0, and thus 


ImQ + ik(3 3 = 0 


(7m + fc 2 /?3/?4) sin 7mft + ik(fo + /? 4 ) 7 m cos 7 m y = 0 

(10) 


(a 2 + k 2 Pi(32) sin a n l + ik{/3\ + fa )a n cos a n z = 0 

( 11 ) 

Substituting for p in equation (1) gives the 
eigenvalues 


kmn — k 7m a n (12) 

The imaginary part of k mn gives the decay rate in 
the x-direction for the mnth mode. Note that 7 m = 
a n = 0 is always a solution to equations (10) and 
(11), but this solution gives the additional conditions 


that 

fc 2 /?3/? 4 + ik(fo + At) = 0 

(13) 

and 

k 2 PiP 2 + ik{Pi +fo) = 0 

(14) 


Thus, for this trivial case, either k = 0 or /3j = /?2 = 
/3 3 = /? 4 = 0. The latter instance represents a plane 
wave solution for a hard-wall duct. 

A general solution of equations (10) through (12) 
is only possible numerically, for example, by using 
a Newton-Raphson scheme (ref. 14). An alternate 
procedure that was used in this study was to use a 
one-dimensional (1-D) finite element formulation of 
a 2-D representation of equation (1) in the x-y plane 
together with boundary condition equations (3) and 
(4) to obtain eigenvalues 7 m . A similar 1-D finite 
element formulation of a 2-D representation of equa- 
tion (1) in the x-z plane together with boundary con- 
dition equations (5) and (6) gives a n . The combined 
eigenvalues k mn are then synthesized, as before, with 
equation (12). This procedure was suggested by 
Watson (ref. 11). 
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Two-Dimensional Finite Element Model 


A 2-D direct solution method for asymmetric 
ducts that need not be uniform in the axial direc- 
tion is described in references 12 and 15. The ad- 
vantage of this formulation is that it is able to model 
the effects of duct liner impedance that vary in the 
axial direction. For the resonator arrays used in 
the present series of tests, this enabled the axially 
varying impedance to be modeled as shown in fig- 
ure 8(d). Here the impedance and equivalent surface 
area of each resonator is modeled separately with the 
equivalent length method, as was done for the modal 
analysis code for the entire resonator array. 

The disadvantage of the particular finite element 
model used in this investigation is that it is necessary 
to assume some dependence of acoustic variables in 
the transverse direction. Usually, for a square-section 
duct such an assumption is not valid unless the so- 
lutions for many cross-modes are summed. In this 
case, however, both side walls are hard and only the 
plane wave mode is cut on, so that the plane wave 
assumption that all acoustic variables are indepen- 
dent of transverse direction is valid in sections of the 
duct where the upper and lower walls are also hard. 
In the duct section containing the resonators, some 
transverse dependence does exist, however. This is 
especially true for the round resonators, since the 
transverse impedance varies because of the resonator 
shape. It was anticipated that the tests described 
in this report would determine the accuracy of these 
simplifying assumptions. 

The finite element scheme employed for this anal- 
ysis was comprised of a matrix of 20 x 150 ele- 
ments. and is fully described in reference 13. There 
the full theoretical scheme for the numerical solution 
for sound propagation in a compressible flow is out- 
lined. For simplicity, this analysis is summarized here 
without the effects of flow. 

The four equations governing the linearized 
acoustic motion in a nonuniform axisymmetric duct 
are 


1 dp 

LOU + — T— - — 0 
p ox 


(15) 


, 1 dp 

LOV + - — = 0 
p or 


(16) 


m 

LOW + — P = 0 
pr 

lop / du dv m v 

~2 + P\ o "X" - d w + - 

c z \ox or r r 


)-' 


(17) 

(18) 


where u.v.w are the acoustic particle velocities, p is 
the acoustic pressure, ~p is the mean density, m is the 
spinning wave number, c is the speed of sound, and 
u> is the source frequency. The cylindrical coordinate 


system of this formulation can be used to model 
a rectangular duct with only a plane wave cut on 
by making m = 0 and r equal to a large number 
(approximating to y in fig. 9(b)), so that w = 0 (i.e., 
plane wave propagating in a large annular duct with 
annulus height much less than duct radius). 

Suppose that within a rectangular element (see 
fig. 9(b)) the acoustic particle velocities u , v, and 
the acoustic pressure p are defined by relations of 
the following type: 


u = (N 1 ,N 2 ,N 3 ,N 4 ){ 


ui 

U2 

U3 

U 4 


(N){u} (19) 


where ui,U 2 ,u 3 ,U 4 are the values of u , v, or p at 
the four node points at the numbered corners of 
the element, and ( N ) are the finite element shape 
functions. 

Define a set of local coordinates (//, f ) within the 
region such that the coordinates of the nodes (%, &) 
are given by (-1,-1), (+1,-1), (+1,+1), (— 1, 4-1) 
for i = l through 4, respectively. The TVj’s of equa- 
tion (19) are shape functions for a linear Serendipity 
element and are defined by 

Ni = + £&)(! + Wi) (20) 

When the simple rectangular element is mapped 
to global (r, x) coordinates it may be distorted to 
fit the required geometry of the problem by the 
transformation 


r={N 1 ,N 2 ,N 3 ,N 4 ){ 


n 

r 2 

r 4 


( N){r } 


x = (N 1 ,N 2 ,N 3 ,N 4 ) 


' X] 
< *2 
Z3 
. X4 


(N){x} (21) 


where (r\,xi), (r 2 , 2 : 2 ) are the global coordinates 
of the corner nodes. The spatial derivative of equa- 
tions (15), (16), (17), and (18) may be expressed in 
local coordinates by 



[ J ]- 1 


J&WI 

U WJ 


where [J] is the Jacobian matrix. 

Using the Galerkin version of the method of 
weighted residuals, we integrate equations (15), (16), 
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(17), and (18) over the element, weighting each equa- 
tion in turn by each of the iVj’s to obtain the matrix 
equations in the nodal variables (ref. 15). It is con- 
venient in integrating these equations to do so in the 
local, rather than in the global, coordinate system 
because this greatly simplifies the process because of 
the convenient limits of integration. Effectively, in- 
tegration is then over a simple rectangular shape of 
figure 9(b). 

Two-Microphone Intensity Measurement 

Technique 

The measurement of intensity at an axial loca- 
tion in a duct can be viewed as an extension of the 
two-microphone technique used to determine acous- 
tic impedances. The analytical basis for the measure- 
ment procedure has been developed by several inves- 
tigators (e.g., see refs. 16, 17, and 18). In essence, 
the method involves a transfer function measurement 
between two points in the standing wave held gen- 
erated by the incident and reflected waves on the 
test specimen. The analysis will be outlined here for 
completeness. 

For plane waves, the acoustic pressures at loca- 
tions x\ and X 2 (see fig. 10) are defined as 


simple relation 


1 + 72 


( 28 ) 


The accuracy of this method depends on the abil- 
ity of the transfer function to provide sufficient in- 
formation to reconstruct the standing wave pattern 
in the tube. The reconstructed wave pattern accu- 
racy depends jointly on the accuracy of the transfer 
function measurement (amplitude ratio and phase 
difference) and on a judicious selection of the two 
microphone locations relative to the test specimen 
surface. Essentially, the method must be capable of 
reconstructing the “true” standing wave pattern in 
the sense that it should be indistinguishable from 
that obtained from a continuous traversing probe 
(see refs. 19 and 20). Clearly, the method relies 
on very accurate amplitude and phase calibration of 
the microphones, which is a tedious and, therefore, 
time-consuming procedure. This constitutes a ma- 
jor disadvantage of the method as reported in refer- 
ences 16, 17, and 18. With this difficulty in mind, 
the two-microphone data were acquired by travers- 
ing a single-probe microphone from one position to 
the other, thereby avoiding the need for tedious 
calibration procedures. 


p( x i) - A (e lkxi + Re lkxi ^j (22) Results and Discussion 


and 

p(x 2 ) = A ( e ikx 2 + Re~ ikx ^ (23) 

Therefore, the transfer function between these two 
locations becomes 

v(xi) fikxi I ikx\ 

P\2 = + (24) 

p(x 2 ) e lkx 2 + Re- lkx i 

Equation (24) can be solved for the reflection factor 
R as 

p ikx\ _ p Jkx 2 

R = 4 ■■ 12 ■■ (25) 

p^~g—ikx 2 g—ikx i x 7 

The acoustic intensity / can then be determined by 

I(xi) = ^Re{p(xi)w*(xi)} (26) 

where v*{xi) is the complex conjugate of the acoustic 
velocity at location aq, as defined by 

v{ Xl ) = — \e lkx ' -Re~ ikxi ] (27) 

pci J 

Finally, the impedance f at the reference plane (look- 
ing from left to right in fig. 10) is determined by the 


Impedance Measurements 

Measurements of normal incidence impedance 
were obtained for both the round resonator arrays 
and the square resonator arrays. Figures 11 and 
12 show the resistance and reactance, respectively, 
versus frequency for a cavity depth of 5.4 in. The 
data shown were taken for resonator 5 in the round 
and square arrays and for an incident sound pres- 
sure level of 120 dB. Data for the round and square 
resonators are depicted by the round and square sym- 
bols, respectively, and have been normalized by pc. 
Figure 12(b) shows segments of the reactance data 
on an expanded scale in the vicinity of the first two 
resonances (i.e. , zero reactance crossings). 

The measured resistances for each of the two res- 
onators are nearly constant with frequency, except 
in the neighborhood of the cavity antiresonances of 
1.26 and 2.52 kHz (see fig. 11). The large increase 
in resistance near cavity antiresonant frequencies is 
characteristic of thin resistive face sheets backed by 
nondissipative cavities. For much smaller diameter 
cavities (i.e., narrow parallel channels) with no re- 
sistive face sheet, acoustic energy dissipation occurs 
mainly in the channel boundary layers. In this case, a 
transmission model (see ref. 21) that accounts for dis- 
tributed dissipation predicts a gradual increase of the 
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channel entrance resistance as the channel antireso- 
nant frequency is approached. In the limit of concen- 
trating all the dissipation in a thin face sheet, a dra- 
matic increase in the effective resistance near the res- 
onator antiresonance is observed, which is consistent 
with the transmission line model. Because the back- 
ing cavities employed in the usual aircraft-type reso- 
nant duct liners have relatively large diameters (and, 
thus, negligible dissipation), a simplified model has 
been employed that adds the impedance of the face 
sheet to the reactance of the backing cavities. This 
model was also employed in the present investigation. 
Therefore, the resistance used in the propagation 
models was approximated by continuous, straight 
lines, faired through the data as shown. The mea- 
sured resistance of the round resonator (which to- 
gether with the face sheet forms a resonator element) 
is about 40 to 50 percent greater than that for the 
square resonator. 

Reactance data over the frequency range 0.5 to 
3.0 kHz are plotted in figure 12(a). As expected, 
the data follow a -cot (kL) trend as indicated by 
the dashed curves in the figure. In the vicinity of 
the antiresonance frequencies, 1.26 and 2.52 kHz. 
-cot (kL) becomes infinite. Resonance frequencies 
are defined by the zero reactance crossings at approx- 
imately 0.61, 1.82, and 3.12 kHz. To better observe 
the small perturbing effect of the fiber metal face 
sheet on the total reactance, the scale of the graph is 
expanded near the resonance frequencies as shown in 
figure 12(b). The fiber metal face sheet contributes 
a small positive reactance ( e lujt time convention), 
as is evident in the figure. By comparing the first 
and third resonant frequencies with the correspond- 
ing —cot (kL) zero crossings, it is evident that this 
effect shifts the resonant frequencies to lower values 
and that the amount of shift increases with frequency. 
These observations are consistent with the small pos- 
itive reactance associated with any thin, porous face 
sheet. 

In view of the 45-percent difference between the 
acoustic resistances of resonator 5 of the round and 
square arrays, dc flow resistance measurements were 
conducted for all seven resonators of the two arrays 
to determine the resistive component variability. As 
described in the “Experimental Setup and Measure- 
ment Procedures” section, these measurements were 
accomplished by attaching a source of controlled air- 
flow to the resonator cavity and measuring the result- 
ing volume flow and pressure drop across the fiber 
metal face sheet. Flow resistance data for all seven 
resonator elements of the round and square arrays 
are shown in figure 13. Note that the smallest in- 
cident velocity was about 70 cm/sec. which corre- 
sponds to a 140-dB incident plane wave. Since the 


main purpose of these measurements was to evalu- 
ate resonator resistance variability, all measured flow 
resistances were extrapolated to zero velocity. 

The zero velocity flow resistances, normalized by 
pc, were plotted versus resonator location for both 
arrays and are presented in figure 14. The dashed 
lines represent average flow resistances of the two 
resonator arrays and serve as references for compar- 
ing the flow resistance variability of the two arrays. 
Note that the average flow resistance for the round 
resonator array is nearly twice that for the square 
resonator array. Also, the variability for the square 
resonator array was much less than that for the round 
resonator array. Visual inspection indicated differ- 
ences in the fiber distribution densities and in the 
wicking of the bonding agent into the active area of 
the fiber metal. Thus, even though the two face 
sheets were taken from the same parent material 
with the intention of having equal resistances for the 
two arrays, significant differences in flow resistances 
were obtained that were apparently due to material 
inhomogeneities and bonding technique. 

Broadband Acoustic Attenuation 

Measurements 

The first series of acoustic attenuation measure- 
ments were performed using broadband noise to de- 
termine the effect of array orientation on the total 
attenuation across the array. Because of measur- 
able differences in the resonator resistances, which 
resulted in an asymmetrical resistance distribution 
with respect to the midspan, or number 4, resonator 
(see fig. 14), attenuation measurements across the 
resonator arrays were conducted for both the for- 
ward and reverse orientations with respect to the 
sound propagation direction. The results are shown 
in figure 15 for the round and square resonator ar- 
rays at uniform cavity depths. Total attenuation was 
obtained by measuring the sound pressure at posi- 
tion B relative to that at position A, as indicated by 
the sketch in the figure. Sound propagates from left 
to right and encounters resonator 1 first in the for- 
ward orientation and resonator 7 first in the reverse 
orientation. Position A is located in the leading edge 
plane of the array and position B is located in the 
trailing edge plane. The attenuation data for each 
resonator orientation were plotted versus frequency 
between 0.3 kHz and 3.0 kHz as shown. The overall 
sound pressure level at the leading edge in all cases 
was 132 dB. 

Figure 15(a) shows the comparison for the round 
resonator array. For this array, there is a small, but 
measurable, change in the attenuation between the 
two orientations, whereas for the square resonator ar- 
ray (fig. 15(b)) the change is barely detectable. These 
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observations are consistent with the greater variabil- 
ity and asymmetry about the midspan resonator in 
the round resonator resistances. A symmetrical vari- 
ability of resonator resistances about the midspan 
resonator would not be expected to produce a change 
in attenuation upon orientation reversal. 

Measured and Calculated Axial Attenuation 

Profiles 

Figure 16 shows (for a representative set of fre- 
quencies) measured and calculated axial attenuation 
profiles over the length of the resonator arrays. These 
data were obtained with discrete frequency sound 
that was incident on the resonator array in the for- 
ward orientation. The cavity depths were set at their 
maximum of 5.4 in. Sound pressure levels for nine 
test frequencies in the neighborhood of the first res- 
onance (610 Hz) were measured at axial distance in- 
tervals of 1 in., starting at the array leading edge. 
As indicated in the figure key, the shaded and open 
symbols represent the calculated and measured pro- 
files, respectively. The calculation results presented 
here were performed with the finite element code, 
so that the complete attenuation profiles, including 
end effects, of the resonator array could be compared 
with the measured profiles. This necessitated that 
the equivalent length impedance distribution model- 
ing procedure, depicted in figure 8(d), be employed 
for the round resonator array. Attenuation rates were 
extracted from the straight-line slopes of the linear 
regions of the profiles as indicated in figure 16. 

For the round resonator array results shown 
in figure 16(a), the agreement between measured 
and calculated attenuation profiles is good to excel- 
lent. with at most a 1-dB underprediction for the 
624-Hz test frequency. The agreement between mea- 
sured and calculated attenuation rates, as deter- 
mined from the linear portion of the profiles (see 
faired lines through data), is excellent. Thus, at least 
for no flow and for impedance nonuniformity separa- 
tion distances up to 0.14 acoustic wavelengths, at- 
tenuation profiles are apparently well predicted by a 
finite element model that allows only axial impedance 
variations. For the square resonator array data pre- 
sented in figure 16(b), the calculated and measured 
attenuation rates are still in good agreement, but 
there are significant differences between the mea- 
sured and calculated profiles. At 500 Hz there is 
a 4-dB difference along the greater part of the array 
length. These differences seem surprising in that the 
measured resistances for the square resonator array 
exhibited much less variability than did those for the 
round resonator array, and the interstitial blockage 
is much less than that for the round resonator array. 
One possible explanation for the greater discrepancy 


is that the square resonator resistances were much 
closer to the optimum value for this duct geometry 
than were the round resonator resistances. It is well 
known that attenuation rates become more sensitive 
to changes in impedance components in the neighbor- 
hood of the optimum impedance values (see ref. 22). 
For this duct geometry, the optimum nondimensional 
impedance at 600 Hz is about 0.16 — 0.14b At this 
frequency, the average nondimensional impedance is 
about 0.88 -I- 0.0« for the round resonator array and 
0.45 + 0.0? for the square resonator array. 

In figure 17, the axial attenuation rates for the 
two resonator arrays are expressed as decibels per 
duct height versus frequency or, alternatively, ver- 
sus the nondimensional parameter kh, where k is 
the wave number and h is the duct height. Thus, 
these results can be applied to any geometrically 
similar duct with the same “equivalent” impedance 
distribution, provided incident energy is carried pre- 
dominantly by the plane wave mode. The measured 
and calculated attenuation rates for the round and 
square resonators are shown by the round and square 
symbols, respectively. The attenuation rates in this 
case were calculated using the infinite waveguide, 
modal analysis, with both the smeared and equiva- 
lent length surface impedance distribution modeling 
procedures (see figs. 8(b) and (c)). The darkened 
symbols represent the smeared impedance distribu- 
tion modeling procedure and the enclosed X symbols 
represent the equivalent length impedance distribu- 
tion model (see fig. 8). The open symbols represent 
the measured attenuation rates as taken from the 
linear regions of the measured attenuation profiles of 
the type presented in figure 16. 

A cursory examination of the data trends reveals 
that the maximum measured and calculated attenu- 
ation rates for the square resonator array are about 
twice those for the round resonator array. This obser- 
vation again confirms that the average resistance of 
the square resonator array was nearer the optimum 
resistance for this particular duct geometry than was 
the round resonator array. Also, the calculated peak 
attenuations (darkened and enclosed X symbols) ap- 
pear shifted to a somewhat higher frequency rela- 
tive to the measured peak attenuations. This fre- 
quency shift of about 50 Hz is most likely due to 
the reactive impedance component being modeled as 
—cot (kL) without including the positive contribution 
from the face sheet. However, the calculated peak at- 
tenuations are within 1 dB of the measured values. 
Also, the two different surface impedance distribu- 
tion modeling procedures result in at most a 1-dB 
difference in attenuation rates, with the equivalent 
length procedure being more conservative. 
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Figure 17(b) shows a similar comparison between 
measured and calculated attenuation rates from the 
finite element analysis method, which accounts for 
impedance variations in the axial direction and. 
therefore, finite length effects. In this calculation, 
the smeared impedance was used to model each indi- 
vidual resonator impedance as an equivalent rectan- 
gular region of impedance, as depicted in figure 8(b). 
As for the equivalent length impedance distribution 
model, this analytical method gives excellent agree- 
ment between measured and calculated results. 

Since the duct height h and the separation dis- 
tance between impedance nonuniformities are nearly 
equal (see fig. 1), the nondimensional frequency pa- 
rameter kh also serves as a nondimensional measure 
of the impedance nonuniformity separation distance. 
The results of figures 16 and 17 suggest that the 
smeared and equivalent length impedance distribu- 
tion modeling procedures are equally reliable for val- 
ues of kh up to at least 0.8. Even though the pre- 
dicted attenuation profiles apparently become less 
accurate as the modeled impedance approaches an 
optimum value, predicted and measured attenuation 
rates are still in reasonably good agreement. This 
circumstance may be due to the increasingly critical 
role of end reflections as the impedance approaches 
optimum. A complicating factor in this argument 
is the greater variability of the round resonator re- 
sistances as compared with those of the square res- 
onator array. However, it is believed that this vari- 
ability can be discounted because very little change in 
the total attenuation spectrum occurs when the array 
orientation is reversed, as shown in figure 15. 

Measured and Calculated Transmission 

Losses 

Next, comparisons of measured and calculated 
transmission losses for three different cavity depth 
configurations of the round resonator array were per- 
formed. No further tests were performed on the 
square resonator array because of inadvertent phys- 
ical damage. Measured transmission losses were ob- 
tained from incident and transmitted net intensi- 
ties as described in the “Experimental Setup and 
Measurement Procedures” and “Analysis” sections. 
The baseline test configuration consisted of all cav- 
ity depths set at a uniform depth of 5.4 in. The 
second configuration consisted of a linear variation 
of the cavity depths varying from 5.4 in. at the lead- 
ing edge to 1.3-in. at the trailing edge. The third 
configuration consisted of the same linear variation 
and orientation, but with the cavity depths reversed 
so that the sound encountered the 1.3-in. cavity first. 
These configurations are depicted in sketches, along 
with the corresponding test results, in figure 18. 


Figure 18(a) shows the measured and calculated 
transmission losses for the uniform cavity depth con- 
figuration. The open and darkened symbols rep- 
resent the measured and calculated values, respec- 
tively, of the transmission losses. The measure- 
ments were obtained at 50-Hz intervals from 0.3 to 
3.0 kHz. Calculations were done at 12 frequencies 
that included the first two resonances. The trends 
for the calculated data follow very well those for the 
measured data. For the second resonance, in the 
neighborhood of 1.8 kHz, the calculated data points 
are shifted about 10 Hz toward higher frequencies 
relative to the measured trend. This difference is 
apparently due to small systematic errors in the 
reactive impedance component. 

For the low frequency, broadband noise suppres- 
sion application envisioned for resonator arrays of 
this type, it is desirable that the passbands associ- 
ated with the antiresonant frequencies be eliminated 
as much as possible. The staggered tuning, or vari- 
able cavity depth, configurations to be discussed next 
are an attempt to accomplish this result. 

Figure 18(b) shows measured and calculated 
transmission losses for a linearly varying cavity depth 
configuration, as depicted by the sketch in the figure. 
Note that less transmission loss variation with fre- 
quency is observed in comparison with the uniform 
depth configuration. Specifically, the transmission 
loss extrema over the frequency range of interest are 
observed to decrease from 1 and 22 dB for the uni- 
form cavity depth configuration to approximately 5 
and 15 dB for the two variable cavity depth configu- 
rations. Thus, the trade-off for increasing the trans- 
mission losses at the minima is apparently a sub- 
stantial reduction at the transmission loss maxima, 
at least for this configuration. 

The calculated transmission loss trend tracks that 
of the measured values, with at most a difference of 
1 dB in the vicinity of the resonant peaks. It is not 
clear why there is more discrepancy between mea- 
sured and calculated results for the variable cavity 
depth configuration than for the uniform depth con- 
figuration. One possibility is that the axial variation 
in the resonator resistances becomes more important 
when an axial variation in the resonator reactances 
is present. 

To further investigate the effect of axial variation 
in the resonator array reactance, the linear variation 
in cavity depths was reversed with respect to the 
sound propagation direction (i.e., the first resonator 
encountered by the sound was set at the minimum 
depth of 1.3 in.). Measured and calculated trans- 
mission losses for this configuration are presented in 
figure 18(c). Again, the general trend of the mea- 
sured transmission losses is well predicted by the 
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theoretical model. However, the calculated transmis- 
sion losses range both above and below the measured 
results by about 1.5 dB. Thus, there is no significant 
difference in the ability of the theoretical model to 
predict the results for the two configurations depicted 
in figures 18(b) and 18(c). The chief difference in the 
overall transmission loss behavior for the two variable 
depth configurations would likely be the substantially 
different reflection experienced by the incident wave 
upon encountering the initial resonator. 

Effects of Mean Flow on Resonator Array 

Attenuation 

The final results to be discussed in this investiga- 
tion concern the effect of mean flow on the round res- 
onator array attenuation. Tests were done at mean 
flow Mach numbers of 0.1 and 0.3 and for the uni- 
form cavity depth configuration. Corresponding flow 
profiles are shown in figure 19. Changes in total at- 
tenuation across the array relative to that for the 
no-flow condition are plotted in figure 20 versus the 
excitation frequency or, alternatively, the nondimen- 
sional frequency parameter kh. Data were obtained 
at 10 discrete test frequencies at approximately equal 
intervals from 0.3 to 3.0 kHz that included both 
resonant and antiresonant frequencies. Attenuation 
changes are referenced to the no-flow attenuation 
as measured with broadband noise (see fig. 15(a)). 
Note that the changes in attenuation are presented 
in terms of the total attenuation across the resonator 
array, as opposed to the attenuation per duct height. 

Clearly, there is a substantial effect of the mean 
flow at the resonant frequencies for the higher flow 
Mach number of 0.3. For the flow Mach number 
of 0.1, the effect is not very significant. It is un- 
likely that these attenuation changes are due to flow- 
induced changes in the resonator impedance since, 
according to reference 23, fiber metal face sheets ex- 
hibit little resistance change for modest mean flows 
of the magnitude examined here. More likely, the 
attenuation change is due to convection effects. 

Of the two propagation models examined in this 
investigation, the finite element model offers more 
flexibility in dealing with realistic geometries. Also, 
the measured and predicted results for the no-flow 
condition were in measurably better agreement than 
for the modal decomposition propagation model. 
Thus additional tests were conducted at the higher 
flow speed ( M = 0.3) to further validate the fi- 
nite element model in the presence of flow. Again, 
these tests were performed with the uniform cavity 
depth configuration. Since the data of figure 16 in- 
dicated a large sound pressure level (SPL) attenua- 
tion across the liner at 624 Hz, that frequency was 
chosen for these comparisons. Figure 21 presents 


analytical predictions for M = 0 and M = ±0.3 in a 
color format that qualitatively depicts the SPL vari- 
ations throughout the test section in the resonator 
array region. As indicated in the legend, the SPL 
increases in intensity as the color changes from the 
black end to the white end of the color spectrum. It 
should also be noted that, for convenience, the res- 
onators are assumed to be on the upper wall in the 
color graphs. The resonator positions are indicated 
to scale by the shaded areas at the top of each color 
contour plot. An intriguing difference between the 
no-flow data and the flow data is the variability of 
SPL across the duct in the presence of flow, as com- 
pared with the constant SPL across the duct at par- 
ticular axial locations with no flow. It is also inter- 
esting to note the small-scale ripples superimposed 
on the larger scale color gradations in the vicinity 
of the resonators. These predicted “hot spots” are 
apparently due to flow effects, as they are entirely 
absent with no flow. Figure 21 also presents quanti- 
tative data in the usual attenuation profile form as 
centerline SPL versus axial distance for each of the 
flow conditions. As would be expected on the basis of 
convection effects, attenuations are less for positive 
flow (i.e., flow with sound) or greater for negative 
flow (i.e., flow against sound), respectively, than for 
no flow. 

To verify that the analytically predicted hot spots 
near the resonator array in the presence of flow were 
real, and not just an artifact of the finite element nu- 
merics, a series of measurements of axial attenuation 
and phase variation were made along a 20-in. length 
of the test section that included the array, and at 
five heights above the array. These results are com- 
pared with the analytical predictions at those par- 
ticular heights in figure 22. Figure 22(a) shows the 
no-flow behavior to provide a baseline for compar- 
isons. Figures 22(b)-(f) shows flow effects in five 
parts, for a flow Mach number of 0.3 and y ranging 
from 0.19 in. to 1.06 in. These data clearly demon- 
strate that the model is capable of predicting not 
only the smooth trends of the no-flow condition, but 
also the spatial variability measured in the presence 
of flow. As seen in figures 22 (b) (f ) , the analytical 
model predicts the experimental measurements in re- 
markable detail for both the SPL and the phase, at 
each of the heights investigated. It should be recalled 
that no changes in the impedance due to flow effects 
are assumed in the propagation model. 

The results presented in this report indicate that 
the modal decomposition model predicts duct liner 
attenuation rates to within about 1 dB per duct 
height for the smeared impedance distribution mod- 
eling procedure for kh values up to at least 0.8. 
The kh value of 0.8 also approximates the largest 
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nondimensional length scale of the axial impedance 
nonuniformities intrinsic to the round resonator ar- 
ray for these particular tests. The resistive com- 
ponent of the round resonator impedance was well 
off the optimum value for plane waves propagat- 
ing in the present duct geometry, whereas the re- 
sistive component for the square resonators was 
nearer the optimum. This circumstance, which 
was unplanned in this experiment, indicates that 
further work should be done in which the round 
resonator impedance is near optimum. Based on lim- 
ited tests with the square resonator, spatial nonuni- 
formities in the impedance may be far more crit- 
ical for the round resonator array near optimum 
impedance values. Comparisons of the finite element 
model predictions with measurements suggest, how- 
ever. that this propagation model should be capa- 
ble of accurately predicting detailed variations in the 
acoustic field near the resonator array in the presence 
of flow and near optimum impedance values. How- 
ever. it was not within the scope of this investigation 
to perform tests at this condition. 

Concluding Remarks 

The objective of this investigation was to eval- 
uate two particular analytical models, together with 
two different impedance distribution modeling proce- 
dures, for predicting the performance at low frequen- 
cies of duct liner configurations built with individu- 
ally tuned resonators. These analytical models had 
previously been developed to predict the high fre- 
quency attenuation of aero-engine inlet duct liners, 
but no satisfactory validation of their accuracy ex- 
isted in the low frequency range. In the low frequency 
application only the plane wave propagation mode 
exists, but the liner impedance is unavoidably charac- 
terized by a spatially varying impedance distribution 
whose length scale may be a significant fraction of the 
acoustic wavelength. This spatial impedance varia- 
tion is a consequence of available structural compo- 
nents and fabrication techniques. 

A '‘smeared’" impedance modeling procedure was 
used to remove the transverse and axial impedance 
variations so that a modal decomposition propaga- 
tion model could be used to predict attenuation rates. 
An equivalent length distribution procedure was used 
to remove the transverse impedance variation so that 
a two-dimensional (2-D) finite element propagation 
model could be used to predict attenuation profiles 
with the axial impedance variations left intact. Some 
of the more significant results of this investigation 
were as follows: 


• Predicted and measured transmission losses for 
frequencies ranging beyond the second resonance 
of the resonator cavity were in very good agree- 
ment for all resonator cavity depths set the same 

• For linearly varying cavity depths, predicted and 
measured transmission losses were still in good 
agreement, but less so than for the uniform cavity 
depths 

• Generally, the transmission loss peaks and valleys 
were less prominent for the array configuration 
with linear cavity depth variation 

• Significant changes in the transmission loss be- 
havior were observed when the array orientation 
(with linearly varying cavity depths) was reversed 
with respect to the sound propagation direction 

Finally, the most interesting and conclusive ev- 
idence for the validity of the finite element model 
was discovered when a set of detailed calculations 
of attenuation and phase profiles near the resonator 
array, with mean flow present, revealed a fine-scale 
structure of acoustic “hot spots” superimposed on 
the general acoustic field structure. Careful mea- 
surements verified the existence of such hot spots, 
which apparently are attributable to convective ef- 
fects alone, since no flow-induced impedance changes 
were assumed in the model. The physics of this phe- 
nomenon is not clear. It is apparently associated with 
flow interaction with the acoustic field by means of 
the nonuniformities in the surface impedance distri- 
bution. Clearly, this is an item for further research. 

The comparison of predicted and measured acous- 
tic performances of the models discussed in this re- 
port shows that both the 2-D modal decomposition 
and 2-D finite element propagation models give good 
agreement with experiment. In addition, it has been 
demonstrated that the finite element model, in par- 
ticular, is capable of predicting, to a remarkable de- 
gree, the detailed variation of sound pressure level 
in the presence of flow along a boundary of varying 
impedance. The work described in this report also 
shows that accurate experimental measurement of 
sound attenuation is a complex task requiring highly 
specialized equipment. Even in a relatively simple 
geometrical configuration, significant errors can be 
introduced when measurements do not include the 
detailed mapping of acoustic flux. 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 
November 17, 1987 
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Figure 2. Sketch of square cavity resonator array. 


17 






QmamL page is 

&nm QUALITY 



Square 
resonator 
array 
installed 
in test 
fixture 


Test 

fixture 


Round 
resonator 
array - 

I 


L-86-255 

(a) Front view. 

Figure 3. Photograph of square and round resonator arrays showing fiber metal face sheet 
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Figure 4. Normal incidence impedance tube for measuring impedances of square resonators. 
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Figure 5. Schematic of normal incidence impedance measurement apparatus and instrumentation. 
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Figure 6. Schematic of experimental arrangement. 
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Figure 7. Acoustic signal processing flow diagram. 
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(d) Equivalent length model (finite element program). 
Figure 8. Impedance models for round resonator array. 
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Figure 11. Measured resistance of round and square cavity resonators (resonator 5 in fig. 6). 
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Figure 12. Measured reactance of round and square cavity resonators. 
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Figure 13. Flow resistance measurements for individual resonator array resistive elements. 
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(b) Square resonator array. 

Figure 15. Attenuation of broadband noise across resonator arrays for two different array orientations. 
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(a) Round resonator array. 

Figure 16. Measured and calculated attenuation profiles (referenced to LE) for selected frequencies. Calcula- 
tions were from finite element analysis method. 






















10 


8 


Smeared 

Equivalent 

Measure- 

Resonator 

impedance 

method 

length 

method 

ment 

geometry 

• 

® 

O 

Round 

■ 

a 

□ 

Square 


6 — 

Attenuation 

rate, 

dB h 

4 - 



a 

□ 




B 


i 


8 

« 


0L 


.3 


.32 


.5 


kh 


.7 


.54 .75 

Frequency, kHz 


.9 


.97 


(a) Modal analysis method. 

Figure 17. Measured and calculated attenuation rates. 
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(b) Linear cavity depth variation, positive slope in propagation direction. 

Figure 18. Continued. 


36 



Transmission 


o 

Measurement 

• 

Finite element 


calculation 

Microphone locations 



for intensity 

Incident 

j- measurements -i 


direction 


ORrPo 


Frequency, kHz 

(c) Linear cavity depth variation, negative slope in propagation direction. 

Figure 18. Concluded. 
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numbers. 
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Figure 20. Change in attenuation because of mean flow for round resonator array. 
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Figure 21. Calculated two-dimensional attenuation contours and centerline axial attenuations, relative to 
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Figure 22. Measured and calculated flow effects on attenuation at various distances from round resonator array 
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